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' Abstract 

In the Frobenius problem we are given a set of coprime, positive integers oi, 02, . . . , Ofe, 
and are interested in the set of positive numbers NR that have no representation by the linear 
form ^jOiXi in nonnegative integers xi,X2, ■ ■ ■ ,Xk- We give a functional relationship that 
completely characterizes the set NR, and apply it to the case when the numbers are in an 
^ , arithmetic progression. 
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Given a set A of coprime, positive integers ai, a2, . . . , a^, one can define the set of positive num- 
bers NR that have no representation by the linear form J2 ^i^i ™ nonnegative integers xi, a;2, • • • , x^. 
It is well-known that NR is finite, and that there is an integer L, such that every integer larger 
r> \ than L is representahle as a nonnegative linear combination of ai, 02, . . . , a^. 

' This is a classic result, attributed to Issai Schur (1875-1941) by Alfred Theodore Brauer (1894- 

1985), who was his doctoral student and teaching assistant at the University of Berlin |^ p. 133]. 
Brauer is also the one who coined the phrase "The Frobenius Problem" as the problem of deter- 
mining the largest element of NR. He named it after Georg Ferdinand Frobenius (1849-1917) who 
would mention it occasionally in his lectures j3- Frobenius, also at the University of Berlin, was 
Schur's thesis advisor. 

Nowadays, the scope of the Frobenius problem is usually extended to include questions relating 
to the entire set of integers that are not representable. Several questions of interest, such as the 
cardinality of the set NR, are related to the so-called Sylvester sums 

Sm — ^ ^ ^ ) 

n&NR 



*An early version of the material in Sections |5| and |3 was published in working paper 11-99 of the Schulich School 
of Business in July 1999. 
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where m is a nonnegative integer. These sums are named in honour of James Joseph Sylvester 
(1814-1897), who posed the question of "How many integers are not representable" for the two- 
variable case [TH] . 

Although the history of the Frobenius problem goes back to at least the late nineteenth century, 
it still is an active area of research, as witnessed by the recent monograph by Ramirez Alfonsm |14j . 
and its bibliography of several hunderd contemporary references. 

In the current paper, we extend the characterization for the two- variable case, derived in |21j . 
to the general case, apply it to the Frobenius problem for arithmetic progressions, and derive 
closed- form expressions for their Sylvester sums. 



2 A Characterization 

For a particular element o of the set A, let TV be the set of positive integers such that, if n is 
representable, then n — a is not. As no two elements of J\f can differ by a multiple of a, and the 
numbers —a -|- 1, . . . , —1 are not representable, it follows that the set TV is a complete, positive 
residue set modulo a, and hence has cardinality |TV| = a — 1. The set TV is a known entity in the 
literature on the Frobenius problem, and first appeared in a paper by Brauer and Shockley |3j. 
Note that NR and TV are equivalent characterizations, in the sense that one uniquely determines 
the other. So, it seems only natural that any characteristic of NR can be described in terms of the 
set TV. This is indeed the case as evidenced by the following theorem. 

Theorem 1 For every function f, the following identity holds: 

E + - /w] = E - E (2-1) 

n€NR neAf n=l 

Proof. Every element of TV, can be represented as nj = j + aruj, j = l,...,a — 1, with rrij a 
nonnegative integer. This gives 

a-1 "ij-l 

^ [f{n + a) - fin)] = ^ E ifi^ + ma + a) - f{j + ma)] , 

n€NR j=l m=0 

where the inner sum is void when mj = 0. The inner sum is a telescoping sum, and simplifies to 

= E if(j + ^j^) - /(j)] = E / w - E /(^)' 

j=l n&N n=l 

completing the proof. □ 



Note that the cardinality of TV is easily derived by taking f{n) as a nonzero constant. 



3 Applications 

With the characterization at our disposal, we are now in a position to derive various properties 
of NR in terms of the set TV. As a first application, take f{n) as a threshold function: /(re) = 1, 
for n > T, and 0, otherwise, setting the threshold r, at the largest element of TV. The right-hand 
side of 1)2. 1|) evaluates to one, so that there must be an re € NR such that n + a > t. It is easily 
seen that increasing the threshold renders the right-hand side of 1)2.11) as zero, so that there cannot 
be an element re € NR, such that n + a > t. This gives Lemma 3 of Brauer and Shockley jSJ: 

L + a = max {n [ re € TV}. (3-1) 
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Theorem 2.3 in Selmer ^7], an expression for the number of integers that are not representable, is 
obtained by taking f{n) = n: 

^o = ^E^-^(«-l)= E K«J' (3-2) 

where the last equaUty follows from the above characterization of the elements of as a residue set. 
The latter equality can also be obtained by taking f{n) = [n/a], where [xj denotes the greatest 
integer that is not larger than x. 

As another application, take f{n) = e"^ and extend the summation on the right-hand side 
of 1)2.11) to include n = 0. Now use the sum formula for a finite geometric series, and divide by 

_ 2^ give the exponential generating formula (egf) of the set NR as 

y e"" = -^—r y e"" (3.3) 

where A/q is the set Af with the zero element added. Note that, by letting z go to infinity, we 
recover ()3.1() . By expanding the exponential in the left-hand side of ()3.3() . and changing the order 
of summation we see that we have also obtained the exponential generating function of the Sylvester 
sums: 

°o m nz 1 

m=0 neNR nSATo J- c i 

In this we can recognize the footprint of the Bernoulli numbers and polynomials through their 
generating function: 

^— = y Bmix)— (3.5) 
e'' — 1 ml 

where the Bernoulli numbers, Bm are simply the Bernoulli polynomials evaluated at zero: Bm = 
Bm{0). Multiplying both sides of ()3.4j) by z, and using (|3.5j) with a change of variables gives an 
expression for the Sylvester sums in terms of the Bernoulli polynomials: 

mSm-i = y Bra{n/a) - B^- (3.6) 

This result can also be obtained directly from Theorem ^ by taking /(n) = Bm{n/a), using the 
difference formula for the Bernoulli polynomials, and Raabe's multiplication theorem. 

4 The Frobenius problem for two variables 

In the classic Frobenius problem, as posed by Sylvester jTSj, the set A consists of two coprime 
integers a and b. This problem has been extensively studied, and a summary of the results can be 
found in |^. Selmer ^| Example 3.1] gives the set Af in the two- variable case: 

J\f= n = 1,... ,a- l} , (4.1) 

and applies 1)3.11) and ()3.2() to give the well-known and classic results L = ab — a — b and 5*0 = 
^(a — l){b — 1). Characterization ()3.6|) gives 

a-l 

mSm-i = a"'-^yB,^inb/a)-B^. (4.2) 

n=0 

While this is an elegant and succinct representation, and determines the Sylvester sums as an 
explicit function of the parameters a and b, it is not the most convenient representation from a 
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computational point of view for large values of a. To arrive at such a representation, one can expand 
the Bernoulli polynomial, change the order of summation, and use the well-known expression for the 
sum-of-powers to derive an explicit expression for Sm- This approach was taken in |^, and further 
details can be found there. Here, we take a different route and take the exponential generating 
function as a starting point. The egf is easily determined using (|3.3j) as 

^f""^. " (e'^' - iVeb^ - 1) ~ 7^^' ^^'^^ 

m=0 ^ ' ^ ' 

Now multiply both sides of this equation by z^, and use the convolution property of the egf,^ with 
the egf of the Bernoulli polynomials to give: 

m(m - l)Sm-2 = E P - mB^^,, (4.4) 

j=0 ^ ^ 

where we have introduced the Bernoullian polynomial (pjix) = Bj[x) — Bj, for notational brevity.^ 
This gives an 0{m?') computational scheme for the Sylvester sums, as opposed to an 0{am) scheme 
using (|4.2)) . To make the dependence upon the parameters a and h more explicit, we can expand 
the Bernoullian polynomials. Since (l)o{x) = 0, we can drop the index j = from the summation, 
replace m by m + 1 , and simplify to give the expression derived by Rodseth ^Hl '■ 

m(m+ l)^jeoV i )\ j ) ' m 

This gives = \{a - - 1), = ^{a - - l)(2a6 - a - 6 - 1), and ^2 = ^ab{a - - 
\){ab — a — b). The reason to repeat these expressions here, in particular the egf (|4.3() . is to showcase 
the simularity with the formulae for the arithmetic case (and its generalization) that we derive in 
the next sections. 



5 The Frobenius problem for an arithmetic progression 

Given the arithmetic progression a, a + d, a + 2d, . . . , a + sd, where a and d are two positive integers 
that are relative prime, one can define the set of positive numbers NR that have no representation 
by the linear form X)f=o('^ + id)xi in nonnegative integers xq, xi,. . . ^Xs- The choice of a = 1 leads 
to the trivial case that the set NR is empty, so we will assume that a > 2. Furthermore, one can 
make the assumption s < a, as otherwise the last element of the arithmetic progression can be 
expressed asa + sd = (ixo + lx(a + (s — a)d), and thus the set NR does not alter if we restrict 
ourselves to the arithmetic progression a, a + d, . . . , o + (s — l)d. Note that, when s < a, there 
are no redundancies in ^ = {a, a + . . . , a + sd}, and none of its elements can be expressed as a 
nonnegative linear combination of the others. 

Roberts showed that the largest element of NR is given by 

+ l^a + (o-l)(d-l), (5.1) 

generalizing the well-known and classic result for the case s = 1, and a Theorem by Brauer for 
the case d = 1. A simplified proof of 1)5. If) was given by Bateman 1_. 

^Convolution property: If the sequence Sm has the egf f(z), and the sequence Tm the egf g{z), then f{z)g(z) is 
the egf of the sequence Sm*Tm = J2T=o (7) 

■^The polynomials 4>j{x) have the egf t^^-fx' '''^'^ ^^^^ actually an older definition of the Bernoulli polynomials, 
that is no longer in common use. Whittaker |22l p. 98] refers to (pii^) ^ the Bernoullian polynomials. Of course, 
apart from a slight notational advantage in the current paper, the difference is immaterial. 



L + 1 
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Grant determined the cardinality of the set NR as 



5n 



a - 2 



+ 1) (a + t) + -(a-l)(d-l), 



(5.2) 



where t = {a — 2) mod s. This generahzes the result by Sylvester [18, for the case s = 1, and the 
result by Nijenhuis and Wilf for d = 1. The proofs of (|5.1|) and (|5.2j) are based upon an explicit 
enumeration of the elements of NR. 

Tripathi [201 determined the set Af as 



A/": 



n — 1 



+ 1 + dn, n = 1, . . . ,a — 1> , 



(5.3) 



and used (|3.1|) and 1)3. 2|1 to derive the above results of Roberts and Grant. 



5.1 The Sylvester sums 

It turns out that it is advantageous to use a slightly different, but equivalent formulation of Af. 
Instead of (|5.3|) . we use 

A/" = |a [n/s] +dn, n = 1, . . . , a - l| . (5.4) 

This is easily seen as an equivalent by virtue of the relationship [x/s] = [(x — l)/sj + 1, for 
all integers x, where [x] denotes the least integer that is not less than x. Combining this with 
characterization (|3.6|) gives: 



a-l 



mSm-i = a™ ^ ^ Bm{\n/s'] + nd/a) - B„ 

n=0 



(5.5) 



Note that for s = 1, we recover formulae (|4.1I) and 1)4.2(1 from the two- variable case. To arrive at 
a computational more convenient formulation, we determine the exponential generating function, 
and start with: 



nGA^O 



a-l 



1 + 



{a\n/s^+nd)z 



a-l 



n=l 

^adz 



1 + 

n=l 

a-l In/s]-! 



ndz 



l + (e" 



\n/s^-l 

1) e'"'^ 

k=0 



odz 



- + (e--l)^ E 



^(ka+nd)z 



n=l k=0 

To evaluate the double sum, we change the order of summation: 



a-l In/s]-! 

E E e^^"''^"''^^^ 

n=l k=0 



a-l 

E E e^'^^'^"-^'^^ = ^ e'^" 

k=0 n=ks+l 
f— 1-1 

I s I \ ^adz 

E 



kaz 



fc=0 



fc=0 
g(fcs+l)d2 



a-l 



Lra=0 



ra=0 



odz 2 e*^-^ 1 



1 e 



ad2 



1 e' 



.dz 



where we have set h = a + sd. Combining all the elements, and simplifying the expression, gives 



E 

nSA/b 



— + (e'^^-l) 



odz 



obz 



1 e- 



(5.6) 
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Now use (|3.4j) to give 

oo 



+ 



m=0 



{e'^'' - l)(e''^ - 1) {e-<^^ - l)(e^^ - 1) 



(5.7) 



Note the simularity to ()4.3() . and, again, observe the footprint of the Bernouhi polynomials. As 
before, in order to invert this expression, we multiply both sides of the equation by z^, and use the 
convolution property of the exponential generating function to give 



0=0 
m 



(5.8) 



Taking account of the fact that the summands are zero for j = 0, and expanding the Bernoullian 
polynomials, gives a generalisation of Rodseth's formula: 



m J 



m{m + l) f^Q^^Q V J + 1 



m + 1\ / j + I 



1 + 



6-1 



1 



m J 



m{m + 1) 



i + l J + l 



(-dr-^-^B„,-jB,.iV 



a - 1 



+ 



-ti+l 



1 



-Bm. 



j=0 i=0 

We note that this gives Sm as a polynomial in \{a — l)/s] of degree m + 2. In particular. 



m 



2So 
125i 

I2S2 



a-1 



+ (2a - 2 + s) 



a-1 



+ (a-l)(d-l), 



-2s{2a + ds) 



a-1 



(6(d - l)a(a - 1) + {3d -2a- ds)s 
s({a + sd){2a + sd) + a'^^ ^ ^ 



+ (ea^ + 3(2a + ds)(s - 1) 
a-1 



a-1 



+ 



+ {a-l){d- l){2ad -a-d-1), 



+ 



a - 1 



+ 



(4a^ + 2(3s - 2)a^ + 2sd{3s - 2)a + 2s^d^{s - 1 
(6{d - l)a^ + 3(2 - s - 2d)a'^ - 3sd{s - 2)a - sd^{s^ - 3s + 1 

(2(2d -l){d- l)a^ - (6^^ - 6d + 2)a^ + 2d{d - s)a + s{l - s)d 
ad{a — l){d — l)[ad — a — d). 



a - 1 



s 

a-1 



+ 



The expressions for Si and S2 in the arithmetic case have not previously appeared in the literature. 
The first expression can be simplified to 5*0 = (a — 1 + r) \{a — l)/s] + (a — l)(d — 1), by repeatedly 
making the substitution s \{a — l)/s\ = a — 1 + s — r, with r = (a — 1) mod s, and is, of course, 
equivalent to (|5.2j) . Unfortunately, employing the same device for Si and ^2 does not yield any 
expression that is much simpler than the ones given. For the arithmetic case, the computationally 
more convenient form of the Sylvester sums does not seem to have an apparent structure nor does it 
exhibit any pleasing pattern, as in the two- variable case. At least, the author has not been able to 
find such. Nevertheless, one should bear in mind that (|5.8|) . although not as aesthetically pleasing 
as (|5.5() . was derived from a purely computational point of view, and for this purpose it is certainly 
adequate. 
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6 The Frobenius problem for generalized arithmetic progressions 



Selmer ^7] generalized the arithmetic progression by taking A = {a, ha + d,ha + 2d, . . . ,ha + sd}, 
with d,h > and gcd(a, d) = 1. As before, we may assume that s < a, as otherwise the last 
element can be expressed as ha + sd = d x a + 1 x (ha + {s — a)d), and thus the set NR does not 
alter if we remove the last element from A. One can also show that, when s < a, there are no 
redundancies in A, and none of its elements can be expressed as a nonnegative linear combination 
of the others. Selmer showed that 



L 





a -2 






s 





and 



-{a-l)[h 



+ d + h-l] + -h 



+ 1 (a — 1 mod s). 



(6.1) 
(6.2) 



Matthews [S| studied the generalized arithmetic sequence in the context of numerical semigroups, 
and determined M as 



A/" 



ha 



n 



+ 11+ dn, n = 1, . . . , a — 1 



(6.3) 



She then applied 1)3. 1|) and 1)3. 2|) to derive Selmer's expressions for L, and the following expression 
for So, that is easily seen to be equivalent to Selmer's: 



So 



1 



-h 



l]{a + t) + ^{a-l){d-l). 



(6.4) 



where t = {a — 2) mod s. 



6.1 The Sylvester sums 

As before, it is advantageous to rewrite (|6.3|) . and take 

M = {ha [n/s] + dn, n = 1, . . . , o — 1} . 



(6.5) 



This allows us, as one would expect from the structure of Af, to leverage the results from the 
previous section. Characterization ()3.6() now gives: 



a-l 



mS, 



m—l 



n=0 



(6.6) 



The computational more convenient form can be derived from the exponential generating function. 
We start with 



a-l 
n=0 



^{ha\n/s~\+nd)z 



a-l 

E' 

n=0 



^{a\n/ s] +nd/h)hz 



and observe that this is the same sum as for the arithmetic case, but with d/h and hz, instead of d 
and z. In the derivation of (|5.6|) we have not made use of the fact that d is an integer, so that we 
can use the result, with the proviso that we replace z by hz, and d by d/h. This gives 



E' 



3(4^1 - 1 



^[ha+sd)z ]^ e~'^'^ - 
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Now use (|3.4j) to give 



^{^ r^i - 1 (e"^" - 1) (e r^i (^'^+^'^)^ - 1) i 



m! (e'^^ - l)(e"^ - 1) (e'*^ - l)(e('^«+^«')^ - l)(e-'^^ - 1) - 1 ' 

To invert this expression, we multiply both sides by z^, and use the convolution property of the 
exponential generating function to give the companion of H6.6|) : 

m(m - l)(m - 2)Sm-z 



m—l / -I 

m — 1 



m 



3 



^ , d^-^-^B^^j_,a^-^<f>j (h [{a - l)/s-] +d) + 



j=0 

E E ( ) ( ^ ) a^-'-'<Pn.-j{h){-dy-'-\ha + sdy-'Bj^,^,{\{a - l)/s] ) 
j=0 i=0 ^ ^ ^ * 



m{m-l)Bm-2- (6.7) 

Note that, for /i = 1, only the index j = m — 1 of the first summation in the double sum remains, 
as 0m(l) is zero for all m, except for m = 1, for which it has the value one. So that, for /i = 1, we 
recover the expression for the arithmetic case. Expanding the Bernoullian polynomials does not 
give any additional insights; we therefore refrain from doing this. 



7 Observations and comments 

We have taken a as an arbitrary element of A. It turns out that this choice, to a large extent, is 
irrelevant. Taking f{n) = x" in ()2.1() gives 

™a _ 1 E ^ ^ „ _ 1 + E ^ ' 

and shows that the left-hand side is an invariant with respect to the actual choice of a. So, if we 
have ai and 02, with corresponding sets 7Vo(ai) and A/'o(a2), then knowledge of one set implies 
knowledge of the other through their generating functions: 

y x^ = —^ y x". 

/^ai Z. — / X^^ 1 — ^ 

Therefore, the most reasonable choice seems to be to take a as the smallest element of A, as this 
minimizes the cardinality of the corresponding set Af. 

The special cases of the Frobenius problem that we considered here have elegant and succinct 
characterizations of the set M. For the general case, however, it is not immediate what its structure 
is. This is probably not an easy question to answer, given that just determining the largest element 
of Af for variable k is NP-hard, as shown by Ramirez Alfonsm JH]. However, for fixed k, polynomial 
time algorithms do exist, see Kannan [S], and algorithms to compute Af have been given by Wilf |23j . 
and Nijenhuis (lOj, among others. 

On a concluding, historical note, in [2T|, we showed that taking /(n) = [n/aj in the two- variable 
case gives 

y[nb/a\ = y l = i(a-l)(6-l), 

n=l neNR 

and provides a new interpretation of this sum as the number of integers that are not representable 
as a nonnegative linear combination of a and b, but omitted to mention that this formula can 
already be found in a paper by Hacks [7j. 
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A Bernoulli numbers and polynomials 

In this section we have collected a number of properties of the Bernoulli numbers and polynomials 
in order not to distract from the flow in the main body where these properties are needed. The 
Bernoulli numbers can be defined by the implicit recurrence relation 

m , _|_ 1 \ 

^p + Ms- =0, for m>l, 

with i?o = 1- This renders the first few Bernoulli numbers as B\ = —1/2, B2 = 1/6, B^ = 0, and 
i?4 = —1/30. Their exponential generating function is given by 



The Bernoulli polynomials can be defined as 



m=U 



m 

m 



r=0 



to give Bq{x) = 1, Bi{x) = x — ^, B2{x) = x"^ — x + ^, B^{x) = x^ — + \x, and B4^{x) 
— as the first fg 

is given by 

E ^™(^)^ = 

m\ — 1 



X — 2x^ + the first few Bernoulli polynomials. Their exponential generating function 



m=0 

The Bernoulli polynomials satisfy the difference formula: Bra{x^\)—Bm{x) = mx™~^, and Raabe's 
multiplication theorem |12j : 

a-l 

Bmix + n/a) = a^~'^Bm(ax). 

n=Q 

Both these are easily proved with the help of the exponential generating function. A more detailed 
expose on the Bernoulli numbers can be found in [3]. 
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